Christensenella minuta interacts with multiple gut bacteria

Introduction Gut microbes form complex networks that significantly influence host health and disease treatment. Interventions with the probiotic bacteria on the gut microbiota have been demonstrated to improve host well-being. As a representative of next-generation probiotics, Christensenella minuta (C. minuta) plays a critical role in regulating energy balance and metabolic homeostasis in human bodies, showing potential in treating metabolic disorders and reducing inflammation. However, interactions of C. minuta with the members of the networked gut microbiota have rarely been explored. Methods In this study, we investigated the impact of C. minuta on fecal microbiota via metagenomic sequencing, focusing on retrieving bacterial strains and coculture assays of C. minuta with associated microbial partners. Results Our results showed that C. minuta intervention significantly reduced the diversity of fecal microorganisms, but specifically enhanced some groups of bacteria, such as Lactobacillaceae. C. minuta selectively enriched bacterial pathways that compensated for its metabolic defects on vitamin B1, B12, serine, and glutamate synthesis. Meanwhile, C. minuta cross-feeds Faecalibacterium prausnitzii and other bacteria via the production of arginine, branched-chain amino acids, fumaric acids and short-chain fatty acids (SCFAs), such as acetic. Both metagenomic data analysis and culture experiments revealed that C. minuta negatively correlated with Klebsiella pneumoniae and 14 other bacterial taxa, while positively correlated with F. prausnitzii. Our results advance our comprehension of C. minuta’s in modulating the gut microbial network. Conclusions C. minuta disrupts the composition of the fecal microbiota. This disturbance is manifested through cross-feeding, nutritional competition, and supplementation of its own metabolic deficiencies, resulting in the specific enrichment or inhibition of the growth of certain bacteria. This study will shed light on the application of C. minuta as a probiotic for effective interventions on gut microbiomes and improvement of host health.


Introduction
The gastrointestinal microbiomes (GMs), as a complex microbial community, are tightly associated with the host health and diseases (Lagier et al., 2018;Dominguez-Bello et al., 2019;de Vos et al., 2022).It is home to hundreds of bacterial species that interact and form networks influenced by diet, physiology, immune response, and probiotics (Hiseni et al., 2021).GMs share the gastrointestinal environment, interact with each other, and associate as a network.Bacterial interactions such as mutualism, commensalism, amensalism or competition (Faust Xu et al. 10.3389/fmicb.2024.1301073Frontiers in Microbiology 02 frontiersin.organd Raes, 2012) are modulated by food intake, host physiological, and immune systems (Weiss et al., 2023), as well as by probiotics intervention (Ma et al., 2020;Zhang et al., 2023).Due to their complexity (Tierney et al., 2019), it becomes challenging to unravel microbial communities' dynamics and functions, particularly in terms of the unknown scope of bacterial interactions (Hibbing et al., 2010).
Metagenomics is a breakthrough technology that allows comprehensive investigation of the genetic material of microbial communities, helping identify and understand distinct species, their functions, and dynamic changes within the communities upon the intricate interplay between environmental factors, interventions, and diseases (Bäckhed et al., 2015;Chu et al., 2021;Shin et al., 2021;Zhang et al., 2021;Ni et al., 2023).Although data-analytical tools such as co-occurrence networking (Mac Aogáin et al., 2021;Narayana et al., 2023) have been used to depict bacterial relations and interactions, the predictions need to be experimentally validated (Jiang et al., 2022).Such validation has been bottlenecked so far by the shortage of cultivated gut microbial resources and also by fastidious workloads (Chang et al., 2019;Liu C. et al., 2021;Abdugheni et al., 2022).The Christensenellaceae family is a bacterial family that has gained increasing attention due to its important role as gut commensal bacteria in human health (Pascal et al., 2017;Lee-Sarwar et al., 2019;Waters and Ley, 2019;Li et al., 2020;Pittayanon et al., 2020;Prochazkova et al., 2021;Tavella et al., 2021;Villanueva-Millan et al., 2022).Early studies revealed that Christensenellaceae was negatively related to body mass index (BMI) and highly heritable in the human guts (Goodrich et al., 2014;Fu et al., 2015;Ruaud et al., 2020;Mazier et al., 2021).Take Christensenella minuta as an example, previous results showed its potential as a probiotic to reduce body weight and inhibit inflammatory reactions (Goodrich et al., 2014;Yang et al., 2018;Mazier et al., 2021;Pan et al., 2022).Thus, C. minuta has significant potential for developing treatment of obesity, metabolic disorders, and inflammatory bowel disease (IBD), such as Crohn's disease and ulcerative colitis (Kropp et al., 2021;Relizani et al., 2022).
A mouse gavage experiment with a C. minuta strain from healthy human gut in some previous research, showed beneficial effects including reducing blood sugar and blood lipids in high-fat dieted mice (Mazier et al., 2021;Pan et al., 2022).In addition, C. minuta produced high levels of acetic acid and moderate levels of butyric acid (Ruaud et al., 2020;Kropp et al., 2021).It has been reported to influence gut microbial diversity and up-regulate the abundance of multiple gut commensal bacteria in previous research, especially for some beneficial bacteria such as Bifidobacterium spp.and Phascolarctobacterium spp.(Mazier et al., 2021;Pan et al., 2022).For instance, it co-occurred with Methanobrevibacter smithii and Oscillospira in multiple population investigations (Goodrich et al., 2014;Klimenko et al., 2018;Ruaud et al., 2020).Similar results were also found in animal experiments: providing C. minuta to germ-free mice transplanted with Christensenellaceae-deficient fecal microbiota from obese humans resulted in the increased abundance of Oscillospira (Goodrich et al., 2014).However, the mechanisms by which C. minuta interacts with other microorganisms remain unclear with few researches focused on cross-feeding via nutrient metabolisms.Ruaud et al. found that C. minuta enhanced the metabolism of M. smithii via H 2 production, which subsequently shifted its fermentation toward acetate instead of butyrate (Ruaud et al., 2020).In addition, C. minuta has limitations in serine, cystine and vitamin B12 biosynthesis (Xiao et al., 2019), suggesting the resources of food nutrition or production by other gut microbes are essential for the growth of C. minuta.In addition, Faecalibacterium prausnitzii and C. minuta have been previously noted to have significant associations in terms of their relative abundances (Otten et al., 2021).However, the specific mechanisms by which C. minuta interacts with other microorganisms remain unclear.
In this study, we enriched the fecal microbiota with or without C. minuta and employed in vitro co-culture with metagenomic sequencing techniques to examine the dynamics of the microbial community in response to C. minuta intervention, and experimentally validated the interaction of C. minuta with predicted bacterial species.The results showed that intervention with C. minuta reduced the microbial diversity but selectively enriched the bacterial metabolic pathways that compensated for its metabolic defects.The positive and negative interactions of C. minuta with other bacteria were mainly governed from nutrient cross-feeding and competition, respectively.Our results provide clues to understand how C. minuta intervention would exert on the gut microbial network and further host health, thus, will shed light on the application of C. minuta as a probiotic for effective interventions on gut microbiomes and improvement of host health.

Culturing and preparation of Christensenella minuta SJ-2
Christensenella minuta strain SJ-2 obtained from the human feces in our previous project was used in this study under an anaerobic environment (Liu C. et al., 2021).SJ-2 was recovered from −80°C in glycerol and cultured in modified liquid Gifu Anaerobic Medium (mmGAM) at 37°C (Liu C. et al., 2021).One and a half mL recovered SJ-2 liquid culture was subsequently inoculated in 100 mL mmGAM and cultured until the optical density (OD) 600 value was 1.0 ± 0.05.This SJ-2 suspension was used for fecal sample enrichment in CMe treatment.

Fecal bacterial enrichment and sample collection
Fecal samples were collected from the healthy volunteer who was confirmed to include no detectable amount of C. minuta by metagenomic sequencing.A 2 g fecal sample was suspended in sterile 0.01 M phosphate buffer solution (PBS) and passed through 70 μm cell sieves (BIOLOGIX, Irvine, United States) to eliminate large particles.Five mL filtrate was considered as "Feces" treatment, and utilized for bacterial isolation and sequencing.One and a half mL of the filtrate was inoculated into either 100 mL mmGAM (GAMe), or SJ-2 suspension prepared in section 2.1 (CMe) in an anaerobic culture bottle, and incubated at 37°C for 10 days.Five bottles were performed for each treatment as five biological replicates.After 10 days of incubation, 5 mL suspension was collected from each bottle and stored at −80°C until analysis.800 μL 0.01 M PBS.Total genomic DNA extraction was performed using the DNeasy PowerSoil Kit following the manufacturer's instructions (Qiagen, Hilden, Germany).A total of 11 metagenomic libraries (including 1 Feace: fecal seed liquid library, 5 GAMe libraries and 5 CMe libraries) were sequenced by Novogene (Tianjin, China) using an Illumina sequencing platform with NovaSeq (PE150, Santa Clara, CA, United States).Each sample was prepared to obtain libraries with a mean insert size of approximately 350 bp.Raw reads were filtered using SOAPnuke v2.1.7 (Chen et al., 2018), and first aligned with the human genome [alignment with Homo sapiens (human)] (RefSeq GCF_000001405.39).The reads aligned to the human genome were removed via bowtie2 v2.3.5.1 (Langmead et al., 2009), and the remaining high-qualified reads in FASTQ format were used for further analysis according to default parameters.
EnrichM v0.6.42 was used to explore the metabolic potential for the MAGs based on Kyoto Encyclopedia of Genes and Genome (KEGG) Orthogroups (KOs) (Kanehisa and Goto, 2000).The non-redundant gene set was constructed using 91 MAGs with software CD-HIT v4.5.8 (Li and Godzik, 2006), resulting in a total of 205,034 non-redundant genes.Subsequently, EnrichM v0.6.4 was employed for annotation, and the differences in community functionality between groups were compared by analyzing gene coverage.The metabolic potential of MAGs was evaluated by annotating modules using EnrichM v0.6.4.MAGs are considered to have complete metabolic capabilities only when the integrity score of Modules is ≥70% (Botté et al., 2023).

Bacterial co-occurrence network analysis
The coverage and relative abundance of all MAGs in each sample were analyzed using coverM v0.6.1 (see footnote 1).Nine total samples were assessed for correlation relatedness as continuous change points using eLsa v1.0.2 (Xia et al., 2011) analysis tools.The computational correlation results were filtered according to Spearman's correlation, and the co-occurrence network was visualized by Cytoscape v3.9.1 (Shannon et al., 2003).

Isolation and identification of gut bacterial species
One mL of fecal filtrate was serially diluted using 0.01 M PBS and 100 μL of the dilution within specified ranges was spread on 150 cm 2 agar plates containing anoxic medium, which were then incubated at 37°C for 21 days.Sixteen different culture conditions were used (Supplementary Table S1).Single colonies appearing on the plates were picked, with 25-30 colonies selected from each plate on the 3rd, 10th, and 21st day, totalling 3 times, and re-streaked onto corresponding agar for purity.The identification of bacterial isolates was achieved via full-length 16S rRNA gene sequencing.Polymerase chain reaction amplicons of the 16S rRNA gene utilizing universal primers 27F (5′-AGAGTTTGATCCTGGCTCAG-3′) and 1492R (5′-GGTTACCTTGTTACGACTT-3′) were sent for Sanger sequencing by Sangon (Qingdao, China) and subsequently identified via EzBioCloud3 (Yoon et al., 2017) and NCBI (National Center for Biotechnology Information 4 ) blast (Federhen, 2012).We used the Neighbor-Joining algorithm with 1,000 bootstrap replicates to construct the phylogenetic evolutionary tree.The multiple sequence alignment of the 16S rRNA gene sequences of 121 isolated strains and one distantly related type strain, Pseudomonas aeruginosa DSM 50071 (NR_026078.1) was performed via the ClustalW function in MEGA11 (Tamura et al., 2021).

In vitro interaction assays
The interactions between the bacterial isolates and SJ-2 were validated using plate confrontation and liquid-culture experiments, respectively, (de Vos et al., 2017;Liu P. et al., 2021).To perform the plate confrontation experiment, 2.5 μL of SJ-2 culture at logarithmic phase was vertically dripped onto a mmGAM agar plate, followed by adding 2.5 μL of another bacterium externally tangentially to the circular droplet.In the experimental setup, both SJ-2 and the target bacterium were subjected to an ecological niche preemption treatment, where each of them was given the opportunity to occupy the niche first (represented as the initial point).We carried out the experiment of dropping the potential interactive bacteria culture vertically first, then dropping SJ-2 at the edge, to verify the influence of one side occupying the ecological niche first on the other side, and the two groups of experiments obtained the same results.The plates were cultivated at 37°C for 2-3 days and photographed to record.Observation and documentation of colony morphology were performed under a stereomicroscope (Nikon corporation SMZ18, Tokyo, Japan) with an eyepiece magnification of 10× and a zoom Xu et al. 10.3389/fmicb.2024.1301073Frontiers in Microbiology 04 frontiersin.orgsetting of 0.75×.For the liquid-culture experiment, SJ-2 cultured in broth medium was centrifuged at 10,000 × g for 15 min after 2 days incubation.The supernatant was filtered using a 0.22 μm filter to remove SJ-2 cells.The conditioned media were prepared by mixing the supernatant with fresh mmGAM medium at equal volume.The tested bacteria were inoculated into the conditioned media at a ratio of 1%.Fresh mmGAM was used to incubate the tested bacteria as a control.The growth rate was determined by measuring the OD at 600 nm under the two conditions during 24 h after inoculation.Due to poor liquid culture conditions and low turbidity of bacterial strains, some bacteria were not subjected to liquid cultivation.Therefore, we performed centrifugal concentration of the bacterial suspension for solid agar-based confrontation experiments.

Mechanism exploration experiment
The physiological and biochemical features of the tested bacteria, including SJ-2 were profiled using ANI MicroPlates (BIOLOG, Hayward, United States) following the manufacturer's instructions.A volume of 3 μL of fresh F. prausnitzii bacterial suspension was inoculated into each well of a 96-well plate (NEST, Tianjin, China) containing 150 μL of preheated culture medium.The plate was then incubated for 24-36 h at 37°C with >95% humidity in a growth curve analyzer (BMG LABTECH SPECTROstar Omega, BMG LABTECH, Ortenberg, Germany) under static incubation conditions and anaerobic environment.The OD 600 was measured approximately every 30 min using a plate reader.Prior to each measurement, the growth curve analyzer was shaken at 700 rpm for 5 min.The growth status of each well culture over time was determined by monitoring the increase in OD 600 (reducing the initial value at 0 h).

Calculation formula for generation time (h). G
Calculation formula for growth rate (/h).

Statistical analysis
The alpha diversity was calculated using the Shannon diversity index, and its significance was assessed by the Mann-Whitney U test.Principal Coordinates Analysis (PCoA) was selected for the analysis of species composition diversity.The Bray-Curtis Index was chosen as the distance method, and PERMANOVA was selected as the statistical method.The intergroup differences in the relative abundance of MAGs, pathways and modules were analyzed for significance using the Linear discriminant analysis Effect Size (Lefse) test, with a significance threshold set at Linear Discriminant Analysis (LDA) ≥ 2.0.β Diversity was used to explore differences between pathway abundance in samples using Principal component analysis (PCA).The potentially interacted bacteria were selected based on the Spearman correlation coefficient between their abundance and the changes in SJ-2 abundance (correlation: SCC ≥ 0.6 or ≤ −0.6).In order to detect changes in the gut microbial ecosystem caused by SJ-2, we removed all the reads assigned as SJ-2 from the dataset and re-normalized them when analyzing community structure and function (Mazier et al., 2021).Mean and standard error of the mean (SEM) were used to analyze and visualize the data.All statistical analysis was performed using IBM ® SPSS ® Statistics 23, including tests for normal distribution and homogeneity of variance.Depending on the data distribution, T-tests, Mann-Whitney's U test, or Kruskal-Wallis test was used to compare the differences and assess the impact of SJ-2 fermentation broth.The effects of substance addiction on F. prausnitzii growth curve were evaluated using repeated one-way ANOVA.To account for multiple testing, p-values were corrected using the Benjamini-Hochberg (BH) method (In all the significance tests for differences, p ≤ 0.05 is considered significant).

Metagenomic sequencing and assembly of gut bacteria
Fecal microbiota was cultivated with mmGAM medium in the presence (labeled as CMe) and absence (labeled as GAMe) of C. minuta SJ-2.The metagenomes from CMe and GAMe, as well as the original seeding fecal samples (labeled as Feces) were extracted and sequenced.Among the total of 11 genomic libraries, 10 were successfully sequenced with one GAMe sample lost.After quality control and filtration, a total of 721,265,291 clean reads were obtained with averages of 72,009,006 (Feces), 77401004.25 ± 2818735.40 (GAMe), and 67930453.60± 1537712.60 (CMe) (Figure 1A, and for the details see Supplementary Table S2).A total of 91 Metagenome-Assembled Genomes (MAGs) were obtained (Supplementary Tables S3, S4), following the criterion of Bowers' classification standard for intermediate to high quality genomes (Bowers et al., 2017).Of the 91 MAGs, 76 were detected in Feces sample, 76 in the GAMe and 65 in the CMe group (Figure 1B, Supplementary Table S3).52 MAGs were shared in all three groups, while 8, 3, and 6 were specifically annotated in Feces, GAMe, and CMe groups, respectively (Figure 1B).In addition, the overlapped MAGs of paired groups were 67, 53, and 58 in Feces-GAMe, Feces-CMe, and CMe-GAMe, respectively.The CMe group exhibits a lower number of MAGs compared to the other two groups.

Christensenella minuta intervention altered the biodiversity and functionality of fecal microbiomes
Both the alpha and beta diversity analyses showed that there was a notable difference between the enriched community and the fecal microbiota.It is not surprising that the growth in mmGAM resulted in changes in fecal microbial diversity, as demonstrated by a decrease in the Shannon index in both the enriched treatments.Additionally, the intervention of SJ-2 further differentiated the fecal microbiomes, with the alpha diversity of the CMe community being significantly lower than that of the GAMe (Supplementary Figure S1).At the family level, Lachnospiraceae was dominant in all 3 groups, while the relative abundance of other families showed apparent changes.The relative abundance of Veillonellaceae, Coriobacteriaceae, and Bifidobacteriaceae increased consistently in both GAMe and CMe groups, while Bacteroidaceae, Selenomonadaceae, and Rikenellaceae decreased in both groups.Meanwhile, Ruminococcaceae and Oscillospiraceae decreased in GAMe but significantly increased in CMe, while Streptococcaceae increased in GAMe but decreased in CMe.Additionally, Clostridiaceae, Atopobiaceae, Lactobacillaceae, Erysipelotrichaceae, and Eubacteriaceae were specifically enriched in the CMe group (Figure 2A).The LefSe results showed that 50 MAGs were specifically enriched in the GAMe samples, while 10 MAGs (excluding SJ-2) specifically represented a high-level redundancy in the CMe group (LDA ≥ 2).Among those MAGs with significant changes, Dialister pneumosintes and Clostridium sp.OF03-18AA were the most representative in the CMe, while Streptococcus pasteurianus, Anaerostipes hadrus and Phocaeicola vulgatus were the most abundant in the GAMe (Figure 2B).
The EnrichM software (Fan et al., 2023) was employed to predict the functional capacity of the fecal microbiota based on the metagenome data (Supplementary Table S5).To gain insights into the KEGG pathway features and differences in the three groups, we performed a functional distribution analysis the results were plotted (Figure 3).Despite the changes in the microbial diversity shown above (Supplementary Figure S1), the overall functionality at pathway level 2 of Feces, GAMe and CMe groups were apparently similar (Figure 3A, Supplementary Table S6).Substance transport, sugar and amino acid metabolism, cofactor and vitamin metabolism, drug resistance, as well as fatty acid synthesis and metabolism were highly enriched in all three groups showing their critical roles during gut microbial survival.These features might suggest that the fecal microbiota possessed functional robustness that enables it to accommodate perturbations in environments.
We further explored metabolic modules in the GAMe and CMe groups.A total of 515 metabolic modules were identified in the GAMe and CMe groups, among which 72 significantly enriched modules were contributed by SJ-2 (Supplementary Table S7).In order to identify the true impact of C. minuta on the fecal microbiome functionalities, the sequences of C. minuta were excluded for further analyses.Performing principal component analysis (PCA) revealed significant differences between the GAMe and CMe groups at the KEGG pathway level 3 (Figure 3B).In general, there were 120 and 101 modules enriched in the GAMe and CMe groups, respectively (Supplementary Table S8).Both the GAMe and CMe groups were abundant in high levels (relative abundances >0.04) of two-component systems, phosphate and amino acid transport systems, phosphotransferase system (PTS), central and other carbohydrate metabolism, ABC-2 type and other transport systems (Supplementary Table S9).The GAMe microbiota particularly enriched PTS and ABC-2 type and other transport systems (Figure 3E, Supplementary Table S9).In the CMe group, the abundances of glutamate, branched-chain amino acids and polar amino acid transport (M00233, M00227, M00223, M00236, M00237), aromatic amino acids (M00040, M00025), and polar amino acid synthesis (M00020) were increased, while polyamine metabolism (M00134, M00133) was enhanced (Figure 3C).Meanwhile, the synthesis of branched-chain amino acids, sulfur-containing amino acids and tryptophan (M00019, M00432, M00570, M00535; M00021, M00035; M00023) showed a significant decrease under SJ-2 intervention (Figure 3C).In the pathways of cofactor and vitamin metabolism, the VB12, VB1, and biotin synthesis modules in the fecal microbiota were significantly enriched under SJ-2 intervention (M00122, M00127, M00577) (Figure 3D).We also observed the enrichment of metabolic pathways such as coenzyme A (M00120), β-oxidation (M00086), and ketone body synthesis (M00088), which regulate the microbiota to a relatively favorable environment (Supplementary Table S10).It seemed that intervention with SJ-2 drove the fecal microbiomes toward the active synthesis of substances that would compensate for SJ-2 metabolic deficiencies such as the synthesis of serine, glutamic acid, tyrosine, biotin, VB12, VB1, thereby enhancing the efficiency of their synthesis in the community after the intervention.On the other hand, the complete metabolic pathways of SJ-2 promoted the  enrichment of microorganisms using arginine, branched-chain amino acids, methionine and tryptophan as metabolic substrates, which were directly supplied through SJ-2 to improve material utilization efficiency (Supplementary Table S11).

Isolation and identification of gut microbes
Totally 1,174 bacterial colonies were isolated from fecal samples using two distinct pretreatment methods and 16 different culture conditions, and further sequenced for their 16S rRNA gene targeting regions greater than 1.4 kb to facilitate their accurate identification.The comprehensive analysis resulted in the classification of 1,174 strains into 121 distinct species, with a stringent threshold of 16S rRNA gene similarity ≥98.7%, which covered 42.6% (29 species) of all the validly named MAGs.However, there were still 22 MAGs that were not explicitly classified into validly published bacterial species.Therefore, we cannot map any isolates to them according to MAG sequences (Figure 5, Supplementary Table S14).Other than the 22 unclassified MAGs, 39 MAGs were not successfully isolated, of which some even showed high abundance in the fecal samples, such as Megamonas funiformis, Roseburia faecis, and F. prausnitzii (relative abundance >3%).In addition, 91 isolated bacterial species were unmapped to qualified MAGs (used for analysis), which accounted for approximately 76% of the total isolates (Figures 5, 6).

Validation of the predicted interactions between Christensenella minuta and other bacterial strains via pair-wise cocultures
In order to validate the predicted correlation in the microbial network between C. minuta and other gut bacterial taxonomies, we tested the interaction of SJ-2 with 16 strains from the bacterial collection of this study, and 7 strains from a previous study (Liu C. et al., 2021) using two methods: (1) pair-wise coculture on agar plates, and (2) sequential cultivation of targeted strain in broth that were conditioned with SJ-2.In the pair-wise coculture, Klebsiella pneumoniae, Phocaeicola vulgatus, Bacteroides caccae, Bifidobacterium bifidum, Bifidobacterium longum, Bacteroides uniformis, Bacteroides fragilis, Bacteroides stercoris, Eubacterium ventriosum, Ruminococcus gnavus, Clostridium innocuum were consistently observed to compete with SJ-2, which was consistent with the prediction in the network [Figure 7A (a-k)].However, Parabacteroides distasonis, Finegoldia magna, Anaerostipes hadrus, Romboussia timonensis, Alistipes putredinis, Flavonifractor plautii, Sellimonas intestinalis, Paraacteroides merdae, Phocaeicola coprocola showed neutral interactions with SJ-2, which cannot correspond their predicted negative relationships in the co-occurrence network ].Meanwhile, among the predicted positive correlated bacteria, Limosillactobacillus fermentum did not exhibit significant growth enhancement when co-cultured with SJ-2; Dorea formigenes competed with SJ-2, which were not consistent with the prediction from the co-occurrence network.Only F. prausnitzii showed relatively increased growth in the presence of SJ-2 [Figure 7A (u-w)].Among all tested strains paired with SJ-2, 12 strains showing consistent tendency with the results in the co-occurrence network prediction indicated the predicted accuracy was about 52.2%.
In comparison, the results of all the liquid culture tests were almost consistent with the relationship predicted in the network except for D. formigenes, including F. plautii, P. distasonis, F. magna, and P. coprocola, which showed neutral relationships in plates but negative correlations in the network (Figure 7B).SJ-2 affected the bacterial growth on either biomass or growth rate or even both.For example, B. fragilis and P. coprocola had prolonged log phases under SJ-2 intervened condition compared with control [Figure 7B (g, t)] while the maximum biomass production of K. pneumoniae, P. vulgatus, P. distasonis, and F. plautii was reduced [Figure 7B (a,b,i,q)].In contrast, both the growth rates and the maximum biomass The co-occurrence network diagram with C. minuta SJ-2 intervention.The pink node labeled with bin.92 represented C. minuta SJ-2, while yellow, green, orange, and purple nodes represented Bacillota, Bacteroidetes, Actinobacteria, and Proteobacteria, respectively.The size of each node corresponded to the status of MAGs in the network, with larger nodes indicating more associations with other bacteria.The connecting lines between nodes were colored red for positive correlations and blue for negative correlations, with thicker lines reflecting stronger correlations.The right side of (Continued) The phylogenetic tree of medium to high quality MAGs was obtained from three different groups.A phylogenetic tree was constructed based on the genome sequences of 91 MAGs using GTDBtk.The tree was built using p_Proteobacteria (GCF_000006765.1)as an outgroup to determine the evolutionary relationships among the MAGs, the outer ring represents the isolation status of the corresponding MAGs at the species level.Cycles from inside to outside: (1) phylum-level assignments (Bacteroidota, orange; Thersmodesulfobacteria, yellow; Proteobacteria, green; Actinobacteriota, pink; Bacillota, blue; Outgroup, gray); (2) isolation/culture status of the 91 (cultivated in this study, purple; strains from previous studies, black; no representative culture, gray).Phylogenetic tree of bacteria isolated and cultured from intestinal samples.Based on the 16S rRNA gene sequences of 121 bacterial strains, a phylogenetic tree (NJ tree) was constructed using MEGA11 with Clustal W algorithm for multiple sequence alignment, using strain Pseudomonas aeruginosa DSM 50071 (NR_026078.1) as an outgroup.The Bar charts were the abundance of cultivated bacterial strains corresponding to MAGs in three treatments (red, Feces; green, GAMe; blue, CMe).Frontiers in Microbiology 13 frontiersin.org3.6 Nutrient competition and metabolic cross-feeding contributed in negative and positive interactions, respectively, between Christensenella minuta and partner bacterial strains In order to explore the interaction mechanism between C. minuta and the other gut bacteria, we performed genome data mining for potential metabolic crossing-feeding and conducted carbon source assimilation tests.Nutrient source assimilation tests showed that there were multiple overlaps in the majority of carbon source utilization between SJ-2 and bacteria showing negative correlation with it in the network, such as B. caccae, B. uniformis, B. fragilis, B. stercoris, P. vulgatus, P. distasonis, B. longum, and B. bifidum resulting in competition between them (Figure 8A).In addition, SJ-2 had a narrower carbon source utilization range than most of the competitors.For example, among all tested 95 nutrient sources, SJ-2 can utilize 25, while B. caccae, B. fragilis, P. distasonis, B. bifidum, and B. longum can utilize over 75% of the substances employed by SJ-2, and at least 20 more that SJ-2 cannot use.
To address the mechanisms of the growth enhancement of F. prausnitzii by SJ-2, we annotated the metabolic pathways of F. prausnitzii and SJ-2 for their genome sequences (Figure 9), and tested ten amino acids as well as fumaric and acetic acids that SJ-2 could apply for their effects on the growth of F. prausnitzii (Figure 8B).The results showed that lysine, proline, tryptophan, cysteine, sodium acetate, and fumaric acid significantly promoted the growth rate, while proline and sodium acetate even increased the biomass production of F. prausnitzii.We observed also that leucine, isoleucine, and arginine prolonged the generation time of F. prausnitzii, with the latter two potentially reducing its overall biomass (Table 1).In addition, high concentrations of methionine, valine, and histidine did not exhibit a significant effect on the growth of F. prausnitzii, highlighting their limited impact.

Discussion
The intestinal microbial community is a complex of beneficial, neutral and potentially pathogenic microorganisms.The intestinal microbial homeostasis was attributed from the interaction of all gut microbes as a network.Christensenella minuta as a potential probiotic increasingly gained attention for its beneficial effects in weight control as well as blood sugar remaining (Kropp et al., 2021;Mazier et al., 2021;Akbuğa-Schön et al., 2023).Based on the current research findings, the colonization of C. minuta has a significant impact on Oscillospira spp., Prevotellaceae, Lactobacillaceae, Erysipelotrichaceae, and Bifidobacteriaceae, which helps in correcting the intestinal microbial displacement caused by obesity (Goodrich et al., 2014;Mazier et al., 2021).However, its interaction with other gut bacteria was rarely explored.In this study, we investigated the interaction of a C. minuta strain SJ-2 with multiple gut bacteria via analyzing metagenomic data collected from in vitro enrichment culture of fecal microbiota, and validated the prediction through one-to-one culturing experiments.After enriched treatment, fecal samples gained noticeable changes in both microbial composition and structure, especially a decrease in biodiversity.The within-sample diversity reflected by alpha diversity was decreased attributing to multiple factors.In addition to the prior microbial disturbance caused by the medium (mmGAM) enrichment, the intervention with SJ-2 further reduced the microbial diversity, particularly for Bacteroidaceae, Selenomonadaceae, and Rikenellaceae.It was not surprising that the presence of SJ-2 consumed some nutrients partially responsible for a discernible decrease in microbial diversity, as evidenced by the high overlap observed between SJ-2 and competitive bacteria such as Bacteroides in the carbon metabolism and substrate utilization investigation.However, some bacterial groups specifically changed with the intervention of SJ-2, such as Lactobacillaceae, Oscillospiraceae, Ruminococcaceae, Atopobiaceae (increased), and Streptococcaceae (decreased), which revealed the specific effects of SJ-2 on fecal microbial community not only realized on the nutrient deficiency.The co-occurrence network analysis results also confirmed that SJ-2 occupied a critical central position that highly associated with most of the gut bacteria, which was consistent with previous studies (Goodrich et al., 2014;Mazier et al., 2021).
Forty MAGs were specifically reduced by SJ-2 in the co-occurrence network, including several opportunistic pathogens such as Finegoldia magna (Shetty et al., 2023), Klebsiella pneumoniae (Lee et al., 2017), Bacteroides fragilis (Goldstein et al., 2017), Bacteroides caccae (Wei et al., 2001;Guo et al., 2016), and Clostridium innocuum (Chia et al., 2017).This was confirmed by functional pathway analysis: the pathogenicity-related genes were reduced in CMe samples.Among all the bacteria that were confirmed to compete with SJ-2, only C. innocuum showed vulnerability in the competition.C. innocuum survives better in lipid-rich environments, and results in creeping fat which is a mysterious feature of Cron's disease (Ha et al., 2020).These combined results reveal that C. minuta's ability for lipid metabolism may contribute to inhibiting the growth of C. innocuum.Furthermore, this result highlights inhibiting the colonization and function of obesity and intestinal disease promoting microorganisms may be one contributor to C. minuta's promotion of weight loss and prevention of intestinal diseases.
In addition, the abundance of the Streptococcaceae family exhibited a specific decrease, with Streptococcus pasteurianus showing the strongest negative correlation with SJ-2 in the co-occurrence network.S. pasteurianus is a potential pathogen for inflammatory diseases, which is known to be enriched in Crohn's disease but reduced as SCFAs accumulated (Nagayama et al., 2020;Wang et al., 2020;Ma et al., 2023).Therefore, the negative correlation between SJ-2 and S. pasteurianus may be attributed to the efficient SCFAs production of C. minuta (Kropp et al., 2021) et al., 2016), which is highly heritable, although accounts for a low proportion of the fecal microbiota (Waters and Ley, 2019;Mazier et al., 2021).Stable colonization has been a challenge for applying most probiotic products.According to the results of this study, a combination of probiotic bacteria such as C. minuta with one or several bacterial strains showing cross-feeding ability may be an option to enhance the colonization and subsequently improve probiotic functions in the human bodies.
Faecalibacterium prausnitzii, which is among the 11 enriched species in the CMe samples, is a potential next-generation probiotic for high butyrate production, anti-inflammation, and prevention of intestinal pathogens (Sokol et al., 2008;Bai et al., 2022).In this study, we demonstrated that SJ-2 cross-feeds F. prausnitzii via producing exogenous acetate, cysteine, proline, and lysine, which were required for the fermentation and reproduction of F. prausnitzii.This phenomenon is consistent with Otten et al's results (Otten et al., 2021), which support the mechanism of the dependency effect of F. prausnitzii on exogenous acetic acid (Duncan et al., 2004).
The intervention of SJ-2 increased the abundance of many potential beneficial bacteria in the family Lactobacillaceae, Oscillospiraceae, Ruminococcaceae, and Atopobiaceae, such as Limosilactobacillus fermentum, F. prausnitzii, and Tractidigestivibacter scatoligenes which strongly revealed their role in host carbohydrate and lipid metabolism regulations.Lactobacillus spp.have been widely used to improve metabolic diseases, enhance gut immunity, and protect gut barrier integrity (Huang et al., 2022).Moreover, members in the Oscillospiraceae and Ruminococcaceae families produce high levels of SCFAs through carbohydrate fermentation, participating in mediating immune response and inducing regulatory T-cell formation (Rooks and Garrett, 2016).In addition, their negative correlation with body weight and inflammatory markers indicated their positive role in host metabolism regulation and antiinflammatory function (Kim et al., 2020;Hu et al., 2021;Leth et al., 2023).Members of Atopobiaceae repress the colonization of drugresistant bacteria in elderly individuals (Ducarmon et al., 2021;Morinaga et al., 2022).Fascinatingly, most of these bacteria contain bile salt hydrolase (BSH) enzymes (EC3.5.1.24)which are involved in antibacterial activity, lipid absorption regulation, and cholesterol reduction processes, akin to C. minuta (Lee et al., 2020;Morinaga et al., 2022;Yang and Wu, 2022;Lin et al., 2023).The positive association between C. minuta and those beneficial bacteria demonstrated that in addition to direct involvement in host metabolism, C. minuta can regulate host's metabolic dysfunctions through collaboration with other beneficial commensal bacteria (Miquel et al., 2013;Leylabadlo et al., 2020;Wu et al., 2021).Demonstration of metabolic cross-feeding between C. minuta SJ-2 and F. prausnitzii.Xu et al. 10.3389/fmicb.2024.1301073Frontiers in Microbiology 17 frontiersin.orgDespite significant changes in the microbial diversity, the fundamental metabolic functions remained stable within the community, related to amino acids, carbohydrates, cofactors and vitamins metabolisms.The intervention of SJ-2 modified the mode of metabolisms of the fecal microbial community, enhancing the transport and synthesis of glutamate, branched-chain amino acid, and polar amino acids, including tyrosine and vitamins B12, and B1.Conversely, it reduced the synthesis pathways of glutathione, methionine, tryptophan, and branched-chain amino acids.Genomic sequence analysis led to the hypothesis that SJ-2 could promote nutritional synergy by cross-feeding with other enriched bacteria.Notably, among SJ-2-enhanced metabolites, biotin is negatively correlated with BMI and is involved in regulating carbohydrate and lipid metabolism, fatty acid synthesis, and amino acid decomposition (Peterson et al., 2020), while thiamin and ascorbic acid are involved in maintaining host immune metabolism, inhibiting oxidative stress-induced NF-κB activation and pro-inflammatory cytokine release, preserving epithelial cell integrity, and reducing host susceptibility (Yadav et al., 2010;Di Renzo et al., 2020).The enrichment of coenzyme A suggests an active β-oxidation process, leading to more ketone bodies and the conversion of acetyl coenzyme A (Canibe et al., 2003).These findings suggest a high level of lipolytic activity in the community, thereby indirectly reflecting the intervention capability (Basolo et al., 2022) and the development potential of C. minuta in lipid metabolism.Moreover, the decreased pathogenicity and dropped abundance of opportunistic pathogens, along with the regulation of other beneficial bacteria, are closely linked to the improvement in the overall health direction of the community.As a result, the microbiota may trend toward a more favorable state for human health.
In summary, we demonstrated the interaction and dynamics between C. minuta SJ-2 and various gut bacteria and tried to explore the metabolic interplay particularly in the cross-feeding between C. minuta SJ-2 and F. prausnitzii.SJ-2 facilitated F. prausnitzii's growth by supplying essential amino acids (lysine, proline, tryptophan, cysteine), acetic acid, and fumaric acid, while reciprocally benefiting from the latter's biosynthesis of nutrients critical for its proliferation (Figure 9).These results enhanced our understanding of C. minuta's probiotic functions and informed strategies for its application.However, this study's limitation to a single sample underscores the necessity for broader investigation involving diverse populations to generalize these interactions.Future work should encompass sampling multiple times throughout the cultivation period to map temporal changes in microbial interactions and explore the influence of varied culture media on these dynamics.In-depth studies at the genetic level will further our understanding of the specific molecular mechanisms underpinning the symbiotic relationships observed.

FIGURE 1
FIGURE 1 DNA sequencing data obtained in this study.(A) Experimental design for resource collections and numbers of metadata and isolates of each treatment.(B) Binning of metagenomic data and assembled MAGs from each sample.

FIGURE 4
FIGURE 4 production of B. caccae, B. longum, B. uniformis, C. innocuum, and F. magna in the conditioned medium were impaired [Figure 7B (c, e, f, k)].On the other hand, both the growth rate and maximum biomass of F. prausnitzii were significantly enhanced with SJ-2 intervention confirming their positive correlation in the network [Figure 7B (w)].

FIGURE 9
FIGURE 9 the network displayed bacteria indirectly associated with C. minuta SJ-2.The cluster at the bottom represented bacteria with no direct correlation with C. minuta SJ-2.The table displayed the correlation of bacteria in the co-occurrence network with C. minuta SJ-2: red, positive; blue, negative.All bacteria in the table were arranged in descending order of correlation intensity.
. This finding offers a new perspective on the development of C. minuta as a probiotic agent for ventriosum; (j) Ruminococcus gnavus; (k) Clostridium innocuum; (l) Parabacteroides distasonis; (m) Finegoldia magna; (n) Anaerostipes hadrus; (o) Romboutsia timonensis; (p) Alistipes putredinis; (q) Flavonifractor plautii; (r) Sellimonas intestinalis; (s) Parabacteroides merdae; (t) Phocaeicola coprocola; (u) Limosilactobacillus fermentum; (v) Dorea formicigenerans; (w) Faecalibacterium prausnitzii (same number as in solid culture medium).regulation of inflammatory bowel diseases.Carbon utilization is overlapped between SJ-2 and various species in the genera Bacteroides and Bifidobacterium, such as Bacteroides uniformis, Bacteroides stercoris, Bifidobacterium longum, and Bifidobacterium bifidum, indicating their negative correlation in the network.Bacteroides spp.and Bifidobacterium spp.are the major commensal bacteria in human gut, most of which show beneficial effects for human health (Zafar and Saier, 2021; Derrien et al., 2022).Although it cannot be concluded that their negative correlation with C. minuta may be harmful to human health, the competition for carbon sources between C. minuta and Bacteroides spp. or Bifidobacterium spp.may result in C. minuta's abundance at a low level in human gut.In fact, previous studies show that B. uniformis is negatively correlated with C. minuta (Upadhyaya the

TABLE 1
Statistical table for maximum growth (absorbance at 600 nm), generation time (h), and growth rate (/h) of F. prausnitzii in culture medium with different doses of additives.